Calculation of the Band Gap Energy and Study of Cross Luminescence in 

Alkaline-Earth Dihalide Crystals. 
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The band gap energy as well as the possibility of cross luminescence processes in alkaline-earth 
dihalide crystals have been calculated using the ab initio Perturbed-Ion (PI) model. The gap is 
calculated in several ways: as a difference between one-electron energy eigenvalues and as a difference 
between total energies of appropriate electronic states of the crystal, both at the HF level and with 
inclusion of Coulomb correlation effects. In order to study the possibility of ocurrence of cross 
luminescence in these materials, the energy difference between the valence band and the outermost 
core band for some representative crystals has been calculated. Both calculated band gap energies 
and cross luminescence predictions compare very well with the available experimental and theoretical 
results. 
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I. INTRODUCTION 

An important revival of the interest in luminescent ma- 
terials is presently observed. This is mainly due to their 
practical applications in several fields, like nuclear spec- 
troscopy, dosymetry, or two-dimensional detectors .used 
for example in medical screens and crystallography^ to 
mention a few. Frequently these luminescent materials 
turn out to be single ionic crystals. Alkaline-earth di- 
halide crystals (pure or doped) are between the most im- 
portant inorganic scintillators. Within this class of mate- 
rials, fluoride crystals with the fluorite structure (CaF2, 
SrF2 and -Bai 7 ^) are the best known scintillators, and re- 
cent work has been devotetLto their study when doping 
with lanthanum®, Mn 2+ BcP and Ce 3+ impurities S~W 
Chlorides are less used materials, though recent theoret- 
ical work on the scintillation properties of SrCl2 ■ Ce 3+ 
samples has been reportedM 1 Experimental studies of thp 
luminescent-emission from impurity centres in SrClj^ 
and BaClJ^ single crystals have also been reported. 
Before formulating a model for the description of doped 
complexes like BaF^ : Ce 3+ , fundamental information 
about the pure crystals should be compiled. Here, as a 
first step towards a more profound comprehension of op- 
tically active impurities, we have focused our interest on 
the description of the pure materials. 

The calculation of the energy band gap (E gap ) of these 
materials is important for several reasons: when doping 
the crystals, it becomes necessary to locate the impurity 
levels in the band gap of the pure crystal. Furthermore, 
the gap is a very important quantity in the first step of 
the scintillation process, namely absorption of radiation 
leading to formation of electron-hole pairs. This step 
influences the global efficiency of the scintillator. Last 
but not the least, it is also important from the theoretical 
point of view, and it gives information on the goodness 
of the model. 

The discovery of the fast luminescent component in the 
emission spectr um , o f BaF^ and the cross-luminescence 
(CL) mechanismli3li3) has brought about a large activity 
in the field of scintillator research.^ BaFi, with a 0.8 
ns CL component is the inorganic scintillator with the 
fastest response. This makes BaF-i a very attractive ma- 



terial for applications in positron emission tomography 
(PET), where a good time resolution is of paramount im- 
portance in order to supress r.ajjidom coincidences and to 
use time-of-flight informationlLj) BaF2 has also been pro- 
posed asA^andidate material for detectors in high-energy 
physicsE3iilP , most recently at the new proton colliders 
SSC and LHC at CERN0 where the important require- 
ment is to differentiate events from different bunch cross- 
ings. In the CL mechanism, an electron is first promoted 
from the core band to the conduction band, leaving a hole 
in the core band, and next a valence band electron recom- 
bines with that hole giving rise to cross-luminescence and 
leaving a hole in the valence band (see Fig. 1). That hole 
can recombine then with the electron promoted to the 
conduction band, leading to selftrapped exciton (STE) 
emission. An important parameter in the study of cross 
luminescence is the energy separation between valence 
and core bands (AEyc)j a $ it determines, for a given 
material, whether the ocurrence of CL is possible or not. 
Specifically, this enenw .difference must be smaller than 
the band gap energyEJTiP (that is, AE C v < E gap ) in 
order to observe CL. This is associated to the fact that, 
if the emission energy is smaller than the band gap, pho- 
ton reabsorption inside the crystal or Auger emission of 
electrons are not possible. 

In the past, there has been a considerable deal of 
work devoted to calculate the band structure of some 
of these pnic, crystals, mainly those with the fluorite 
structure.EJO'' The purpose of the research carried out 
here is to study pure AX2 crystals, where A stands for 
Mg, Ca, Sr, Ba and X for F, CI, Br, I, by using in- 
stead a cluster approach. |3jo this end we have used the 
Perturbed Ion (PI) model.O) The present work is an ex- 
tension of our own previous research on ionic crystals, 
where first we calculated the band gap energies of alkali 
halide crystals] 2 ^ and then studied several properties of 
the scintillators NaI:Tl+ and CsI:Tl+.EP Given the great 
technological significance of cross luminescence processes, 
we have also calculated the energy difference between the 
valence band and the upmost core band for some AX2 
materials, and compared it with the relevant gap energy. 
This allows to predict the occurrence of cross lumines- 
cence in these solids. Cluster approaches have been used 
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in the past to study the CL mechanism in alkaline-earth 
fluorides by Andriessen et a/.EP , and also by Ikeda et 
alB& 

The paper is structured as follows: In the next section 
we show how to apply the PI model to the calculation 
of gap energies between valence and conduction bands 
and between valence and core bands. The results of the 
calculations are presented and discussed in section III. 
Section IV presents our conclusions. 

II. CALCULATIONAL METHOD 
A. Brief summary of the theoretical model 

A detailed description of the PI method has been given 
in our previous work on band gap energies of alkali halide 
crystals^ and in the original worksE3'E2P , so here we just 
give a brief resume. In order to calculate the electronic 
structure of ionic solids, and in particular the band gap 
energy, the traditional band theory approach could be di- 
rectly applied. For doped crystals, however, the transla- 
tional symmetry is lost, and the Bloch's theorem does not 
apply. Cluster models, which have been successfully ap- 
plied to the analysis of local crystal properties, provide an 
alternative, although any cluster approach has some in- 
herent difficulties in dealing with delocalized conduction 
states (see, for example, ref. |34|). The cluster approxima- 
tion is based on a partition of the solid into a "cluster 
region" and "the rest of the crystal" . The main issues one 
has to address when using this approximation are the fol- 
lowing : a) the election of the best partition of the crystal 
into "cluster" and "host" regions; b) a proper quantum 
mechanical description of the finite cluster; c) a precise 
description of the ions surrounding the cluster (this is 
necessary for a proper account of cluster-lattice interac- 
tions); and d) the consistency between the descriptions 
of the finite cluster and of the surroundings. As a matter 
of fact, issues c and d have received less attention than 
issues a and b. Issue d is frequently neglected, and the 
"rest of the crystal" is often simulated just by using point 
charges.© More accurate lattice models have beenCjOn- 
sidered, for example by Winter, Pitzer and Temple ,£30' 
who introduced effective core potentials to describe the 
lattice cations nearest to the cluster. The significance of 
the clusterJattice consistency was stressed by Kunz and 
coworkers Ji!fP who included the so-called localizing po- 
tentials in the description of the cluster. Besides, issue 
a is also a delicate point, as it is difficult to avoid sur- 
face effects at the cluster boundary when studying bulk 
properties with cluster models. 

In the PI model, the "cluster" is reduced to its min- 
imum size, a single ion, and cluster-lattice interactions 
are described in the framework of the The ory o f Elec- 
tronic Separability (TES) of Huzinaga et a/JM'EJ Being 
the "cluster" a single ion, boundary effects are avoided 
right from the start. Furthermore, the cluster approxi- 
mation can be rigourously formulated within the TES, 
as cluster-lattice orthogonality is a fundamental require- 
ment of that theory. Each ion in the lattice acts on the 
"cluster" density through an effective potential which in- 



cludes nuclear, Hartree and exchange contributions. The 
electronic structure of each inequivalent ion in the crystal 
is self-consistently determined, thus avoiding the use of 
empirical parameters. EZP Another feature of the PI model 
is that it does not invoke the Molecular Orbital LCAO 
approximation. The one-center character of the model 
leads to a large computational-time saving compared to 
any multi-center cluster approach. 

The atomic-like orbitals used to describe each ion are 
expanded— into a large set of Slater type (STO) basis 
functionsEilB* because of their superior performance. 
The outputs of the PI model are a set of crystal-adapted 
wave functions for each ion, which are fully consistent, 
in the Self-Consistent-Field (SCF) sense, with the lattice 
potential, and the effective energies of the corresponding 
ions in the field of the crystal lattice. The effective energy 
is written as a sum of the net energies (or self-energies) 
of the ions in the crystal and the ion-lattice interaction 
energies. Once those quantities have been obtained, the 
total crystal energy is obtained as a sum of monocentric 
terms, involving net energies, and bicentric-cpntributions 
coming from the interaction between ions.ElP A correla- 
tion energy correction following the unrelaxed Coulomb- 
Hartree-Fock- (uCHF) prescriptioncija) is added to the 
net energies.© The PI method is particularly well suited 
to model an impurity center in a crystal because we can 
describe the ions surrounding the impurity with lattice- 
consistent wave functions, rather than with Hartree- Fock 
free-ion wave functions. Applications of the PI model to 
these problems can be consulted in refs. |2S|^ . 

B. Band gaps 

The band gap energy can be defeaed as the difference 
between the two following energieseJ) a) the energy nec- 
essary to ionize the crystal by removing an electron from 
the top of the valence band (the threshold energy E t ); b) 
the energy gained by putting an electron at the bottom 
of the empty conduction band of the perfect crystal (the 
electronic affinity x)'- 

E gap =E t - X - (1) 

To calculate x we have to model an electron in the con- 
duction band, that is in a delocalized state. The strong 
ion-lattice orthogonality required by the TES would force 
that electron to be localized on a given ion. We can, 
however, obtain a reasonable description of a delocalized 
conduction state by relaxing some of the orthogonality 
requirements imposed by the TES (see below). The elec- 
tronic affinity x is always a smalL quantity, close to zero 
for all these ionic materials JlSila) and our description 
gives values of the correct order of magnitude (a few 
tenths of eV). As typical errors in measured gaps are ~ 
0.5 eV, the quantitative determination of the band gap 
energy is not critically affected by our approximations. 
All the calculations, ate performed at the experimental 
crystal structures-LElL^r 

We proceed now to calculate the band gap energy in 
several ways. In a first calculation we approximate the 
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threshold energy E t by the energy eigenvalue (with op- 
posite sign) of the highest occupied anionic orbital ob- 
tained in a PI calculation for the pure crystal at the HF 
level. This corresponds to the Koopmans' approximation 
to the ionization potential. To obtain x> we simulate a 
singly charged alkaline-earth cation A + as an impurity in 
the field created by the pure crystal (A + : A 2+ X 2 ), and 
equate x t° the eigenvalue (with opposite sign) of the 
outermost s-orbital of that cation. To allow for the delo- 
calization of the orbital over a substantial region of the 
crystaJ-we-,enlarge the basis set taken from Clementi and 
Roette!rt3 ) for A + with some diffuse functions and relax 
the orthogonality requirement between that s orbital and 
the surrounding lattice. This eigenvalue is close to zero 
(in no case is — e greater than 0.2 eV). According to the 
definition (eq. 1), our first calculation of the gap reads: 

Egap — £ ns (A + : A 2+ X 2 ) — e mp (X : A 2+ X 2 ) (2) 

where n=3,4,5,6 (m=2,3,4,5) for Mg, Ca, Sr and Ba 
(F, CI, Br and /L,|Cespectively. In recent studies, de 
Boer and de GrootoH 1 have shown that the conduction 
band of those ionic crystals formed with alkali or alkaline- 
earth cations is associated mainly to anionic levels. Thus, 
we have checked the effect of centering the conduction 
electron on an anionic site. As long as we enlarge the 
anion basis set with enough diffuse functions, the value 
of x is almost insensitive to the election of specific center 
for the conduction electron. 

In a second calculation, the band gap energy is ob- 
tained as a difference of total energies, as in a typical 
ASCF calculation. This approach includes orbital relax- 
ation effects, improving thus over Koopmans' approxima- 
tion. The independent processes of removing an electron 
from a halogen anion and of placing that electron on the 
conduction band read now: 

Et = E crysta i(X° : A 2+ X 2 ) — E crysta i (A 2 +X~), (3) 

and 

X = E crystal (A 2+ X 2 -) - E crystal (A+ : A 2+ X 2 ), (4) 

respectively. E crysta i(X° : A 2+ X 2 ) represents the en- 
ergy of the crystal with a single neutral halogen impurity 
and E cryst ai (A + : A 2+ X 2 ) the energy of a crystal with 
a delocalized electron although centered on an alkaline- 
earth ion (see above). In the process of calculating the 
energies of the systems X° : A 2+ X 2 and A + : A 2+ X 2 , 
the ions surrounding the impurities are described with 
the lattice-consistent wave functions obtained in the PI 
calculation of the pure crystal. With these assumptions, 
E t and x reduce to a difference of effective energies 

E t = E eff (X° : A 2 +X 2 ) - E eff (X- : A 2 +X 2 ) (5) 

X = E eff (A 2 + : A 2 +X 2 ) E eff (A+ : A 2+ X 2 ). (6) 

E e ff(X° : A 2+ X 2 ) is the effective energy of a neutral 
halogen atom in an otherwise perfect crystal A 2+ X 2: 
and a similar interpretation holds for the other effective 
energies (for details, consult refs. ^and ^6|). E t gives the 



main contribution to E gap , while x only provides a small 
correction of magnitude ~ 0.2 eV, as in the previous 
calculation. The ASCF calculations are performed at the 
HF and uCHF levels; the second ones include, besides 
orbital relaxation effects, correlation corrections to the 
HF energy. 

C. Cross luminescence 

We extend now the previous expressions to the calcu- 
lation of AEvc, the energy difference between valence 
and core bands. In some AXi crystals (see section III.C 
below), the core band nearest to the valence band is as- 
sociated to the outermost p-orbitals of the A 2+ cation, 
but in some others it is associated to the outermost s- 
orbitals of the X~~ anion (when the anionic s-eigenvalue 
is less negative than the cationic p-eigenvalue) . Our first 
calculation is again a simple difference between energy 
eigenvalues (Koopmans' approximation). For the first 
case: 

AE vc = e mp (X-)-e n _ hp (A 2+ ), (7) 

and in the second case: 

AEyc = £m P {X~) — e ms (X ), (8) 

with a notation consistent with eq.(2). We can improve 
the calculation of AEvc allowing for the relaxation of 
the hole states created in the curse of the CL process 
(see Fig. 1). We obtain first the core binding energy, 
Ec- Allowing for orbital relaxation, the expression for 
Ec is similar to eq.(5): 

E C = E eff {A* + : A 2+ X 2 ) E eff (A 2+ : A^X 2 ), (9) 

if the core band is formed by the outermost cationic p- 
orbitals, or 

E c = £ e// (X°(ms W) : A 2+ X 2 ) - E eff (X~ : A 2+ X 2 ), 

(10) 

if the core band is formed by the outermost anionic s- 
orbitals (the notation X a (ms 1 mp 6 ) indicates the elec- 
tronic configuration of the halogen species after removing 
one electron from the ms orbital) . Ec is the energy neces- 
sary to remove an electron from the core band out of the 
crystal. By substracting from this energy the threshold 
energy of the crystal, given in eq.(5), we obtain the fol- 
lowing approximation for the energy difference between 
core and valence bands: 

AE VC = E C -E t . (11) 

As before, the ASCF calculations are carried out at the 
HF and uCHF levels of theory. 

III. RESULTS AND DISCUSSION 

An important part of the whole process involves to 
solve the electronic structure of the pure crystal. The 
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PI model is a general method for dealing with crys- 
talline compounds of any spatial group, and the only 
inputs required are the spatial group and the lattice con- 
stants. CaF2, SVF2, BaF2 and SrCli adopt the fluorite 
structure (Fm3m in the international notation for spa- 
tial groups), MgF2 a rutile-type structure (P^/mnm), 
CaCl-2 and CaBr2 adopt distorted rutile structures 
(Pnnm), MgCl-2 has the cadmium chloride structure 
(i?3m), and MgBr2, Mgl^ and Ca/2 the cadmium io- 
dide structure (P3ml). Finally BaCl 2 , SrBr2,rJ3aBr 2 
and Bal 2 adopt the P6C7 2 structure (Pbnm)!^ The 
treatment of complicated crystalline structures is not a 
challenge for the PI calculations. Experimental geome- 
tries have been obtained from ref. |5(], except those of 
BaC'h, BaBr2 and Pa/2, which have been taken from 
refs. |^, |5^ and |5^, respectively. In the Pbnm structures 
the two anions are in slightly nonequivalent positions, so 
they have different eigenvalues and effective energies. To 
calculate the threshold energy in these cases we will take 
out from the crystal the least bound p-electron, which 
corresponds to an anion X~ in a definite site. In all the 
other cases, anions are in equivalent positions. 

A. Band gaps 

Results of the calculated band gaps of fifteen alkaline- 
earth dihalide crystals are presented in Table I. We give 
the calculated gaps at three diferent levels of theory and 
compare these with experimental results whenever they 
are available. We have collecled— the experimental re- 
sults from several sources. Ej'OlZIIIj'LEP The AenF gaps 
from eq. (2) show the expected behavior in ionic solids, 
namely, the band gap decreases when moving down the 
periodic table along the alkaline-earth (halogen) column, 
leaving fixed the halogen (alkaline-earth) ion. The in- 
fluence of a different spatial symmetry does not affect 
at all these trends, so it seems that its effect on band 
gap energies is small. Comparison to experimental val- 
ues shows a systematic overestimation of the gap. The 
gaps obtained by the ASCF(HF) procedure are smaller 
than those calculated by substracting eigenvalues. Re- 
laxation of the halogen atom in response to the removal 
of one electron from the anion (which is included in E L 
at this level but not in the previous calculation) is the 
main reason responsible for this effect. The influence 
of orbital relaxation on x is rather negligible. We ob- 
serve for the ASCF(HF) method an underestimation of 
the band gap energies of fluorides but still an overesti- 
mation for chlorides. Inclusion of Coulomb correlation 
effects is crucial to achieve quantitatively good results, 
as the ASCF(uCHF) calculations show. Upon inclusion 
of correlation, the band gap energies change in the cor- 
rect direction both for fluorides and chlorides, increasing 
the gap in the first group and decreasing it in the sec- 
ond group. The larger response of the anionic density 
cloud to correlation effects is also responsible for this im- 
provement. On the other hand, \ remains practically 
unchanged, as before. 

Theoretical studies piX2aF2 and MgF2 have been car- 
ried out by Catti et alsBlEy using CRYSTAL, an ab initio 



periodic Hartree-Fock program. They have obtained 
ground state properties, like equilibrium geometries or 
elastic constants, and also have shown the calculated 
band structure for these crystals. The band gap ener- 
gies obtained are 21 eV and 20 eV for CaP 2 and MgF2, 
respectively. The HF approximation tends to overesti- 
mate the gap in ionic crystals E3 At first sight, it seems 
that these results should be compared to our Ashf cal- 
culation, which, although still overestimating the gap, 
gives closer agreement with experiment. However, both 
calculations should not be directly compared. In refs. 

the gap is identified with the energy difference be- 
tween the highest occupied and the lowest unnoccupied 
orbitals, that is a LUM O — HOMO calculation. In our 
Achf calculation, while E t can be identified still with 
the HOMO, x ca n not be identified with the LUMO 
anymore, as it is the eigenvalue (with opposite sign) of 
an occupied orbital. We think that this fact, together 
with the appropriate inclusion of diffuse basis functions 
in the description of the conduction state to simulate de- 
localization (diffuse valence orbitals on cations were not 
included in the above mentioned works) are the main rea- 
sons for our improvement over their results. This shows 
that even an energy difference between eigenvalues can 
be a quite good approximation when the conduction state 
is more realistically described- 
More recently, Ikeda et alx3 have calculated the band 
gaps of CaF2, SrF2 and Bap2 by using active clusters of 
different sizes embedded in a field of point charges rep- 
resenting the Coulomb potential created by the crystals. 
The electronic structure of those active clusters is solved 
by the use of the discrete variational (DV) Xa method.Eif 
The band gaps obtained by Ikeda et al. show some de- 
pendence on the size of the active cluster employed, but 
for a fixed cluster size, show the correct trend, that is a 
systematic decrease of the band gap when changing the 
cation size from Ca 2+ to Ba 2+ . A meaningful quanti- 
tative comparison with our set of results is not direct, 
because their calculation identifies the gap with a differ- 
ence between one-particle electronic levels (like our Aehf 
calculation), but including correlation effects at the same 
time footing as exchange through the Xa model. Nev- 
ertheless, the gaps they obtain with the (FyA4) + cluster 
are almost identical to our most accurate ASCF(uCHF) 
results. 



B. Cross luminescence 

Now we turn to the calculation of the energy differ- 
ences involved in cross luminescence processes. We show 
in Table II the energy eigenvalues (with opposite sign) 
of the core orbitals nearest to the valence band, for 
four representative crystals. We note that the outermost 
core band is formed by cationic alkaline-earth p-states in 
SVF2, BaF2 and BaCl2, and by anionic halide s-states 
in SrCh- So the energy of the core band level (Pc) is, 
in the Koopmans' approximation, identified with the 3s 
eigenvalue of CP in SrCfo, with the 5p eigenvalue of 
Ba 2+ in BaF2 and BaCl2 and with the 4p eigenvalue of 
Sr 2+ in SrF2- The energy differences AEvc are shown 
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at three levels of theory in Table III, together with exper- 
imental results for fluorides as given in ref. H, The Ae HF 
calculation underestimates AEyc m these materials, be- 
cause the Koopmans' approximation overestimates the 
threshold energy E t more than Ec- We can appreciate 
this item more clearly when passing from the Aehf to 
the ASCF(HF) calculation. Here, contrary to the band 
gap calculation (see Table I), AEcv increases in SrF 2 , 
BaF 2 and BaCl 2 , showing that anionic orbital relaxation 
effects are more important than cationic ones: both E t 
and Ec are lower than in the AenF calculation, but the 
lowering of E t is more important, resulting in an increas- 
ing of AEcv- In SrCl 2 , such behaviour is not observed, 
and this is due to the fact that the core band is an anionic 
band. The ASCF(HF) result lowers slightly the magni- 
tude of AEyc for this crystal. This effect is very small, 
however, showing that the error in applying Koopmans' 
approximation is nearly the same for Et and Ec in this 
case. The ASCF(HF) results include orbital relaxation, 
and the different relaxations of core and valence holes are 
responsible for all the effects discussed. Finally, the effect 
of Coulomb correlation (ASCF(uCHF) column) is, ex- 
cept for BaCl 2 , very small, showing that the corrections 
introduced in the calculation of E v and Ec are similar, 
and tend to cancel out in the energy difference. Then it 
seems that correlation tends to shift occupied bands by 
a rather constant amount with respect to the conduction 
band, leaving energy differences between occupied bands 
unaffected, at least for fluorides. 

In the fifth column of Table III we show the possibility 
of occurrence of CL in these crystals. As we explained 
in the introduction, if the energy difference AEyc be- 
tween the valence and the outermost core band is smaller 
than the energy gap E gap between conduction and va- 
lence bands, then it is possible to observe CL. This con- 
dition holds in BaF 2 , so we can assert that,CL is possible, 
as it has been experimentally observed. t3t3 Regarding 
SrF 2 , we obtain that AEyc is larger than E gap by 1.7 
eV, so cross luminescence is not predicted by our cal- 
culations. The situation is even more unfavourable for 
SrCl 2 and BaCl 2 crystals (for which we obtain a value 
for AEvc larger than the corresponding band gap en- 
ergy by 8 eV) as well as for the rest of AX 2 crystals not 
shown in table III. The combination of fluorine and bar- 
ium, which gives the most ionic compound between those 
studied in this work, gives a favourable limiting case for 
CL applications. Barium is the cation with the smallest 
ionization potential in the alkaline-earth series, so it is 
easier to remove an electron from the 5p orbital of Ba 2+ 
than from other alkaline-earth cations. This leads to a 
lower value of Ec for barium compounds. On the other 
hand, fluorides show the largest band gap energies be- 
tween alkaline-£arth dihalides systems (as in simple alkali 
halide crystalsEfP ) , because it is more difficult to remove 
an electron from a F~ anion than from any other halide 
anion. As a consequence, the difference in eq.(ll) is the 
lowest in magnitude. 

The cluster model calculations performed by Ikeda et 
al.cr lead to the same conclusions, although the specific 
value obtained for AEyc again depends a little on the 
size of the active cluster. They obtain a band gap for 



BaF2 which is 3 eV larger than AEyc, so their calcu- 
lations predict also the occurrence of cross luminescence 
for BaF2 . The band gap of CaF2 is smaller than the en- 
ergy difference between the valence and core bands. For 
SrF2, the values of E gap and AEyc are similar, again in 
agreement with our results. 

IV. CONCLUSIONS 

We have calculated the band gap energies of fifteen 
alkaline-earth dihalide AX 2 crystals. For some of these 
crystals we have also calculated the energy difference be- 
tween core and valence bands in order to study the pos- 
sibility of cross luminescence (CL) processes. We have 
followed two different ways to estimate these quantities: 
first, we have taken differences between energy eigenval- 
ues, and next we have carried out ASCF calculations, 
both at the HF level and with inclusion of correlation 
effects. To this end we have used the Perturbed Ion (PI) 
model, supplemented with some additional assumptions 
to deal with an electron at the bottom of the conduction 
band. This is a cluster-type model which achieves full 
lattice-cluster consistency. Correlation effectsJiave been 
included using a model proposed by ClementiJifP Within 
the PI model we can deal with crystals of any spatial 
point group symmetry without any problem. The signif- 
icance of an accurate estimate of the band gap energy in 
these crystals as a preliminary requirement for the study 
of luminescence properties has been pointed out, as well 
as the importance of the energy difference between core 
and valence bands in the CL process. The calculation of 
the last quantity is possible because our calculation is an 
all-electron calculation and does not invoke any approx- 
imations like frozen-core treatments. Orbital relaxation 
of the core hole is also allowed. However, one has to rec- 
ognize that our method gives no dispersion for the core 
or valence bands so the applicability is restricted to cases 
when these bands are narrow enough. 

Overall we have obtained a rather good estimation of 
the energy band gaps, with a computational effort which 
is much less than that required by any other cluster multi- 
center approach or by standard band structure methods. 
Gaps compare quantitatively well with experimental re- 
sults whenever they are available, and favourably to pre- 
vious theoretical calculations for fluorides. The CL mech- 
anism has been predicted as possible only in BaF 2 , in 
agreement with the experimental observations. 
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Captions of figures 
Figure 1 

Schematic view of the cross luminescence process. 
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TABLE I. Calculated band gaps, compared to available experimental values. SPG gives the spatial point group of the 
crystal, in the international notation. AenF is a difference between one-particle orbital energy eigenvalues. ASCF refers to a 
difference between the total energies of the crystal in appropriate electronic states (see text). All energies are given in eV. 





Crystal 


SPG 




ASCI (HI ) 


ASCI (uCHl ) 


Exp. 




MgF 2 


P4: 2 /mnm 


14.3 


10.8 


12.8 


12.4° 




CaF 2 


Fm3m 


12.7 


10.0 


11.9 


12. 6 




SrF 2 


Fm3m 


12.4 


9.7 


11.5 


11. 1 6 




BaF 2 


Fm3m 


12.0 


9.2 


11.0 


10. 5 6 




MgCh 


R3m 


10.7 


9.2 


7.9 


7.5 C 




CaCl 2 


Pnnm 


10.2 


9.1 


8.0 


6.9 C 




SrCh 


Fm3m 


9.6 


8.5 


7.4 


7.5 C 




BaCh 


Pbnra 


8.9 


7.8 


6.8 


7.0 C 




MgBr 2 


P3ml 


10.5 


9.3 


7.6 






CaBr 2 


Pnnm 


9.7 


9.0 


6.7 






SrBr 2 


Pbnm 


8.5 


8.0 


6.5 






BaBr 2 


Pnma 


8.3 


7.7 


6.0 






Mgh 


P3ml 


9.9 


9.2 


7.2 






Cal 2 


P3ml 


9.2 


8.4 


6.3 






Bal 2 


Pbnm 


8.1 


7.4 


5.4 





a : Ref. 57 
6 : Ref. 46 
c : Ref. 58 



TABLE II. One-particle orbital energy eigenvalues for the outermost p-orbital of the A 2+ cation and the outermost s-orbital 
of the X~ anion. All energies are given in eV. 



Crystal 




e(X-) 


SrF 2 


23.6 (4p) 


36.9 (2s) 


BaF 2 


17.7 (5p) 


36.3 (2s) 


SrCh 


27.0 (4p) 


25.6 (3s) 


BaCh 


20.4 (5p) 


24.9 (3s) 



TABLE III. Energy differences between the upmost core band level and the valence band level calculated from different 
theoretical methods (as in Table I), and possibility of cross luminescence (CL). All energies in eV. 





Crystal 


AtHF 


ASCF(HF) 


ASCF(uCHF) 


CL 


Exp. 






SrF 2 


11.2 


13.2 


13.2 


no 


? 






BaF 2 


5.7 


7.8 


7.7 


yes 


yes 




SrCh 


16.0 


15.8 


15.4 


no 








BaCh 


11.4 


12.2 


14.8 


no 
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